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LOCAL FACTORISATION FOR THE DYNAMICS OF QUANTUM SPIN 

SYSTEMS 

SVEN BACHMANN AND ANDREAS BLUHM 


Abstract. Motivated by the study of area laws for the entanglement entropy of gapped ground 
states of quantum spin systems and their stability, we prove that the unitary cocycle gener¬ 
ated by a local time-dependent Hamiltonian can be approximated, for any finite set A, by a 
tensor product of the corresponding unitaries in X and its complement, multiplied by a dy¬ 
namics strictly supported in the neighbourhood of the surface dX. The error decays almost 
exponentially in the size of the neighbourhood and grows with the square of the area 


1. Introduction 

The general locality result for the dynamics of quantum lattice systems proved in the seminal 
work of Lieb and Robinson [T] has been refined and extended in many directions, and it has 
found countless applications. The Lieb-Robinson bound is best understood as expressing the 
linear growth of the support of local observables, and the Lieb-Robinson velocity is therefore 
referred to as the ‘speed of sound’ of the quantum lattice system. In this letter, we use it to prove 
a strong version of locality, namely directly at the level of the unitary cocycle implementing 
the dynamics in any finite volume, the result being uniform in the volume. The Hamiltonian 
generating the dynamics can be time-dependent. Indeed, the result applies not only to the 
actual time evolution of the system, but, in fact more importantly, to the spectral flow studied 
in [2], and used in [3]. The latter article motivates the present work, which provides a simple 
proof of Lemma 1 therein. 

Concretely, let A be a finite subset of the lattice and A C A. Let be the unitary 
evolution generated by the family We prove that can be approximated in norm by 

the product of two unitaries: firstly the purely factorised dynamics <8> secondly a 

‘surface operator’ s)* supported on the R-fattening of the surface dX of X. The error 

decays faster than any inverse power in R and grows like |5Ap. Beyond its structural interest, 
the result and its particular dependence on the surface of X plays an important role in the 
understanding of the stability of area laws for the entanglement entropy, see again [3] as well 
as [1]. The factorisation presented here, which holds in any dimension, should also be put in 
parallel with the similar factorisation of the ground state projections proved in mi. We note 
in particular that the error there grows similarly with the square of the surface of X. 


2. Setting 


We consider a quantum spin system defined on a countable set F, which is equipped with its 
graph distance d{x,y) > 0 for any G F. We assume that there is a F : [0,oo) —>■ (0,oo), 
satisfying the two conditions 


( 1 ) 

( 2 ) 


||F|| := sup'^ F{d{x,y)) < oo. 


xer 


y&r 


snp ^ < oo. 


F{d{x,z)) 


Let ^ : [0, oo) —>■ (0, oo) be a non-increasing function that is logarithmically superadditive, 
namely 


( 3 ) 


C{a + b) > C{a)C{b). 
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Then the function F^{r) := F{r)^{r) satisfies ([2|) again, as 

i{d{x,y))i{d{y,z)) ^ ^{d{x,y)^{d{y, z) ^ ^ 

i{d{x,z)) ~ i{d{x,y) + d{y,z)) ~ 

for a constant < Cp- We assume moreover that 

lim ^(r)r" —>■ 0 
r—>-oo 

for all n G N, in which case F^ satisfies ©• Furthermore, the function 

C(i?):= e(r), i?GNU{0}, 

r>R+l 


is strictly decreasing and decays faster than any inverse power. Indeed, for any n € N, 


lim C{R)R'^ < lim 
R—^oo R—^oo 


i sup 

r>ij+l 




(r - I)*" 


= 0 . 


To any finite A C T, we associate the finite dimensional Hilbert space R/^ = ®x£aT~Lx, where 
diniHx < oo. We denote = F{R\) the set of observables supported in A. It is canonically 
embedded in .Aa', A' D A by tensoring with the identity. 

Let T € (0, oo) and —T < s,t <T. The Hamiltonian F[\{t) € A.a is a sum of time-dependent 
local interactions, 

ZcA 

with T(Z, t) = T(Z, t)* € Az- We assume that the map [—T, T] 3 t T(Z, t) is continuous 
for all finite Z. Furthermore, we assume that 


'J'lk := 


1 

sup -^ 

x,yer F^{d{x,y)) 


E 

Z3x,y 


sup 

t&[-T,T] 


||T(Z,t)|| < OO. 


The self-adjoint F[\{t) generates the unitary cocycle U\{t,s) G A.a solution of 

I^LaCL s) = HAit)UA{t, s), Ua{s, s) = 1. 

The associated Heisenberg evolution 

:=UA{t,srAUA{t,s), 

for any A G A.a, is such that = rjf {i[HA{t), A]) and that 

(4) 4’*orr(A) = ri’^(A). 


The dynamics satisfies a Lieb-Robinson bound of the form 


( 5 ) 




2PIIPII 

< " " min 

Cf 


- s) ^^d^ddix,y))] 

xGX yGY 


where A G Ax, B G Ay and 


did) = 


ex' 


p(ug|r|)-l ifd(X,y)>0 


exp(uj|r|) otherwise 


see [ 7 ], and [2] for the time-dependent case. The positive constant = 2||T||^C^ is the Lieb- 
Robinson velocity. 

Let X C F be a finite subset of F. The inner boundary of X is defined by 


dX:={xeX-. d{x,V\X) = 1} 
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For any i? G N U {0}, consider the sets 

:= {x G r : d{x,X) < R} 

:= {x G r : d{x, T\X)<R} 

OrX :=X^nx^, 

and note that dX = \ Xn . 

We shall finally assume that there is G > 0 such that for any finite X C F and all i?, 
(6) \dRX\ < GR’^mi 

where | • | denotes the volume of the set, i.e. its cardinality. 


3. The factorisation 

In this section, we consider X a fixed finite set such that the map 

N 9 n i-A |X„ \ 2£.n-i\ ^ 0 

is non-increasing. For any i? G N U {0} and A any finite set such that X C A, let 

Mr := {Y C diR/ 2 ]X : y n X / 0 and y n A \ X ^ 0}, 

and we shall write R/2 instead of for simplicity in the following. With this, we define a 

surface energy by 

S{R,t):= T(y,t). 

ZgMr 

Note that 

(7) sup ||S(ii,t)||< ^ sup ||T(Z,t)||<|a^/ 2 ^|||Fg||||T||^ 

x,y&dR/2X ZeMfl;x,yGZ*e[-T,T] 

<G||F^||||T||5(i?/2)'='|5X|, 

by ([6]). Finally, let UQ^x{i-,s) G Aq^^x be the differentiable unitary cocycle defined as the 
solution of 

= Tg’^xiSi^A))Uaj,x{t,s), Ua,ixis,s) = 1 . 

_^ 

Theorem 1. For a finite X C F as above, i? G N, and any finite A C F such that X C A, 


C4(i, s) - {Ux{t, s) (g) U/^\x{t, s)) UaRxit, s)* 


<|5X|||T||^||F|||t-s 


2C{R/2) + \dX\K {R/2f^C (-R/2) 


where k 


G^IIFj 


Proof. By unitarity, 

UA{t, s) - {Ux{t, s) g UA\x{t, s)) Uanx{t, s)* 


V{t,sr -ila^xiPsY , 


where 

V{t, s) = UaY, s)* {Ux{t, s) g UA\x{t, s)) . 

We note that V{s, s) = 1 and that 

s) = -UaY, sYHAit) {Uxit, s) g UA\xit, s)) 

+ UaY, sY{Hx{t) + H^\x{t)) {Ux{t, s) g U^\x{t, s)) 
= T^/{K{t))V{t,s) 


1 ) 
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where 


m = - 


E 




ZCA: 

20X^0 and ZnA\X^0 

By the standard estimate 

V{t,s)* -Ua^x{t,sr = V{t,srUa^x{t,s)-l 

d 


< 


dr 


V{r,s)*Uajix(r,s) 


dr 


< |t - s| sup TX’^(it'(r)) + Tg;^^^(S’(i?,r)) 

rG[t,d 

it suffices to control the sum of the generators, uniformly in r and A. 

Claim 1. For any i? € N: 

sup \\K{r) + SiR,r)\\<2\dXmkmmm-, 

r£[s,t] 

Claim 2. For any ii € N, and for any A € 




Cp 


' A 

We first prove Claim 1. By definition, 

X(r) + 5(ii,r) = - 

the norm of which can be estimated as 


< ^^\\A\\\dX\{R/2fi{R/2) (e^«l"-^l -1). 


E 


ZcA:ZnA\9H/2X^0, 

znA:^0,znA\A:^0 


snp \\K{r) + S{R,r)\\<Y^ ^ EE sup ||^'(Z,r)|| 


re[s,t] 


nGN xex„\x^_, L 


+ x(n>i?/2 + l) E E sup l|'h(Z,r)|| 


Multiplying and dividing by F^{d{x, z)), we indeed obtain 

sup ||A:(r) + S{R, r)|| < mkWn E \ + xin > R/2 + l)e(n)) 

„6N 

<2\dxmk\\naR/2). 

In the last inequality, we used that, by assumption on X, 

\x^\x,-i\<\x,\x,\ = \dx\. 

We now prove Claim 2, which is a well-known consequence of the Lieb-Robinson bound in 
the case of a time independent Hamiltonian. Let A G Ay with Y C A' C A. First, we note that 

= -i Ha{s),t/{A) 

This follows immediately from the cocycle property ([3]), 

° = r^A {i[HAk),r/{A)]) + rf (^) 

and the fact that is an automorphism. With this, 
d 


U,S 


oT^f (A) = tI’"' (i[LfAN - HA'{u),rA^iA)]) , 
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and since — T^f{A) = JJ ° d?’, we obtain 

-'rA/^(^)|| < [ ||[^(^>«)>'^a’^(^)]|| dn 

zeA\A' ZBz'^^ 

<^;^-(e"‘"'-''-l) E E “P lltlZ.ujll^J^FjMi,!,)) 

^ ^ 2eA\A'Z32“^b-d x£Zy£Y 

by the Lieb-Robinson bound ([5]). In order to get an upper bound, we replace '^z^zYlxez by 
SxgaSzsojz! multiply and divide by F^{d{x,z)) and use the convolution condition to carry 
out the sum over x to finally get 


||rP*(A)-TjA-4)|| J2F^{d{y,z)) 

^ 2 eA\A' y€Y 

< ^|S|FU(d(y,A\A'))(e’^d--d _ 1). 

Claim 2 follows by setting Y = djij 2 X and A' = OrX, so that 

|y|eKE,A\A')) < \dx\G{R/2faRm 


by ([6]). 

Finally, since 




Claim 1 and Claim 2 jointly with ([ 7 ]) yield the statement. 


□ 


One may be interested in a similar factorisation with the pure tensor product on the right of 
a surface unitary, namely of the form Udj^xit, s)* [Ux{t,s) 0 Ca\A'(C'S))- The analog of V{t,s) 
being here W{t, s) = [Ux{t, s) 0 Ca\a:(C •s)) UA{t, s)*, the necessary cancellations of interaction 
terms occurs only if one considers the adjoint evolution, namely taking the derivative with 
respect to s. We have 

-i^W{t,s) = W{t,s)a^jCiK{s)), 

where cj^^(A) = U\{t, s)AUA{t, s)*. In that case, let tjQ^x{t,s)* be the solution of 

-'^^Uaax{t,s)* = a^f^^{S{R,s))Ua^x{t,s)\ Ua^xitA)* = 1 , 

so that 


( 8 ) UA{t,s)-Ua^x{t,sr{Ux{t,s)0UA\x{t,s)) = W{t,sr-Ud^xit,sy 


< 


aY{K{r))+aY^^iSiR,r)) 


t^r 


dr. 


From there, the proof runs along similar lines as above. Claim 1 remains unchanged. As for 
Claim 2, we note that (T^*(A) — (A) = — JJ o cryf{A) dtt so that 

||o-a'*(^)-^A*(^)|| ^ Z] / ll[^(^’“)’^A'*(^)]|| dw. 

2ga\A' zbz"^ ^ 


It remains to note that || [T(Z, u), (A)] || = || u)). A] || and observe that the Lieb- 

Robinson bound ([5]) is symmetric in its arguments to conclude that the upper bound of Theo¬ 
rem [1] holds for ([8|) as well. 
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